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The  cross-correlation  technique  for  signal  detection  and  parameter  estimation  is  widely  used  in  communica¬ 
tion  and  radar  systems.  For  popular  binary  waveforms,  the  properties  of  their  ambiguity  functions  have  been 
studied  for  years.  However,  the  results  are  generally  restricted  to  specific  codes  which  are  either  deterministic  or 
pseudorandom.  From  the  uncooperative  reception  viewpoint,  without  a  priori  knowledge  of  the  emitted  signal  in 
many  cases,  the  cross-correlating  signals  may  be  considered  truly  random.  The  report  investigates  the  ambiguity 
function  of  truly  random  binary-phase-coded  waveforms  as  an  approximation  to  those  waveforms  commonly 
employed  in  binary-modulated  pseudorandom  systems/encoded  radar  systems. 

Here,  in  a  statistical  sense,  the  ambiguity  function  of  a  truly  random  binary-phase-coded  waveform  is 
analytically  derived  in  which  the  normally  used  deterministic  cross-correlation  process  is  replaced  by  its  ensemble 
average.  Various  doppler  responses  are  presented  and  discussed.  The  results  are  compared  with  those  obtained 
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by  transmitting  an  aperiodic  maximum-length  pseudorandom  sequence.  It  is  shown  that  the  ambiguity  function  of 
the  latter  case  is  closely  represented  by  the  ensemble-average  response  of  the  truly  random  binary  signal  of  equal 
length.  With  the  transmitted  signal  unknown  to  a  cross-correlator,  the  characteristics  of  the  ambiguity  function 
derived  here  provides  predictive  characteristics  useful  for  practical  cross-correlator  system  design. 
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THE  AMBIGUITY  FUNCTION  OF 
RANDOM  BINARY-PHASE-CODED  WAVEFORMS 


INTRODUCTION 

Since  the  early  development  of  simple  continuous  wave  (CW)  and  pulse  radar  systems,  the  radar 
waveform  concept  associated  with  matched  filter,  pulse-compression,  or  coded-waveform  techniques  has 
been  one  of  the  most  critical  elements  in  radar  design  (11.  It  is  well  known  [1,21  that  the  use  of  the 
matched  filter  in  radar  systems  can  be  considered  as  optimum  predetection  processing  for  white  Gaus¬ 
sian  noise  interference,  and  the  output  of  the  filter  is  indeed  the  noise-corrupted  input  signal  cross- 
correlated  with  a  replica  of  the  transmitted  signal.  The  matched  filter  receiver  is  mathematically 
equivalent  to  a  cross-correlation  receiver,  or  simply  a  cross-correlator.  Generally,  the  detected  signal 
out  of  the  matched  filter  or  the  cross-correlator  is  characterized  by  the  so-called  "ambiguity  function’ 
which  is  determined  by  the  transmitted  waveform. 

The  cross-correlation  technique  for  signal  detection  and  parameter  estimation  is  widely  used.  A 
typical  example  is  the  detection  of  the  emitted  signal  presented  at  two  sensors  remotely  located  in  a 
transmission  medium  [3],  Interesting  problems  may  include  the  geolocation  of  the  transmitter  by 
cross-correlating  the  signals  received  over  two  separate  paths.  In  this  case,  the  extracted  parameters  of 
interest  would  be  time  difference  of  arrival  (TDOA)  and  differential  doppler  (DD).  In  general,  precise 
extraction  of  parameters  requires  the  received  signals  to  be  cross-correlated  over  a  two-dimensional 
ambiguity  surface  on  TDOA  and  DD  domains.  Obviously,  predictive  measures  of  the  observer's  reso¬ 
lution,  accuracy,  processing  gain,  and  interference  rejection  will  be  provided  if  the  ambiguity  function 
of  the  received  signal  can  be  determined. 

For  popular  discrete-coded  waveforms,  the  properties  of  their  ambiguity  functions  have  been 
studied  for  years.  The  results  are  generally  restricted  to  specific  codes  either  deterministic  or  pseu¬ 
dorandom  such  as  Barker  codes,  maximum-length  shift-register  sequences,  and  others  [1,4],  However, 
in  many  cases  without  a  priori  knowledge  of  the  transmitted  signal,  one  would  like  to  have  similar 
results  for  cross-correlating  waveforms  incurred  in  the  binary-modulated  pseudonoise  systems  or 
encoded  radar  systems. 

Indeed,  in  order  to  simultaneously  produce  a  high  data  link  and  reduce  the  bit  energy,  a  longer 
transmission  is  practically  applied.  From  the  reception  viewpoint,  the  discretely  coded  waveform  of 
finite  code  length  can  thus  be  considered  aperiodic.  Here  we  considered  random  binary-phase-coded 
waveforms  as  an  approximation  to  these  pseudorandom  waveforms  and  replaced  the  cross-correlation 
processing  by  its  ensemble  average. 

In  this  report,  the  ambiguity  function  of  a  random  binary-phase-coded  signal,  which  is  considered 
aperiodic  and  finite,  is  statistically  formulated  and  discussed.  In  addition,  the  responses  of  various 
doppler  Alters  are  obtained  and  compared  to  those  generated  through  an  aperiodic  maximum-length 
pseudorandom  sequence  (5). 
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FORMULATION  OF  THE  AMBIGUITY  FUNCTION  OF 
RANDOM  BINARY-PHASE-CODED  WAVEFORMS 

In  radar  or  digital  communication  systems,  the  transmitted  signal  in  complex  representation  is 
generally  given  by  sit,f0)  —  uit)  exp  (J2w  f0t )  where  /0  is  the  carrier  frequency  and  uit)  is  the  com¬ 
plex  transmitted  modulation  with  uit )  —  ait)  exp  (Jmit)).  Here  ait)  and  mil)  are  the  amplitude  and 
phase/frequency  modulation.  For  discrete  binary-phase-coded  waveforms  of  interest,  uit)  can  then  be 

expressed  by 


uit)  -  P„  (r)  exp  0  8,).0<  r  <  N  ft 


where  AT  is  a  finite  integer,  -  0  or  w,6  is  the  pulsewidth,  and  P„it)  -  1  with  (n-1)  8  <  t  <  n  8  and 
*  —  1, 2, ...  ,N.  For  convenience,  the  expression  exp  (J  $„)  is  denoted  by  c„  throughout  the  report. 
Clearly,  {cn}  indicates  a  sequence  of  l’s  and  -l’s.  Nt  and  1/8  are  the  code  length  and  the  pulse 
repetition  frequency. 

When  the  doppler  effect  can  be  treated  as  a  simple  shift  of  the  carrier,  it  is  easy  to  obtain  the 
ambiguity  function  through  the  cross-correlation  of  with  sit*  r  ,fof),  where  r  and  /  are  the 

time  delay  and  the  doppler  shift  Generally,  the  function  is  characterized  by 

X  ir,f)  —  J_m  u'it)uit*  t)  exp  (—  J2  n  ft)  dt  (2) 

where  u'it)  is  the  complex  conjugate  of  uit).  For  the  binary-phase-coded  waveform,  the  ambiguity 
function  can  be  obtained  by  directly  substituting  Eq.  (1)  into  the  above  equation.  The  result  is  [1]: 


Xir ,/)  -  Xf  (f ,/)  £  c„cm  +  k  exp  (-  J2  w  fin  -  1)8) 

R-l 

N-(k+ 1) 

+  X,  («-■?,/)-  £  c„c,  +  *  +  iexp(-;2ir/n8)  (3) 

i?»i 

where  -  00  </<  °°,  r  —  Ac  8  +  f  with  0  <  r  <  8  and  *  —  0, 1, 2, ...  N  -  1,  and 

Xp  it.f)  -  J*Q  exp  (—  jin  ft)  dt. 

The  generalized  form  of  Eq.  (3)  is  of  particular  interest  since  it  includes  both  continuous-  and 
discrete-time  cases  and  can  be  applied  to  an  aperiodic  waveform.  As  cited  in  the  first  section  of  this 
report,  we  consider  {c„}  as  an  approximation  to  pseudorandom  sequences,  truly  random,  and  redefine 
the  ambiguity  function  from  a  statistical  viewpoint.  In  other  words,  we  assume 

£  |x  c*j  “  &  |£  CI  fy|  "  •  •  “  o,  U  t*j)  and  E  c,2  -  E  £  c,4  -  l.and  replace 

Xir  ,f)  of  Eq.  (2)  by  its  ensemble  average,  i.e.,  X2ir,f)  —  E  [X  (t  ,/)  X'irJ)].  By  directly  substitut¬ 
ing  the  result  of  Eq.  (3)  into  X2ir,f)  and  replacing  jin  f  8  by  y,  one  obtains 


N-k  N-k 


chch  *  +  l^m  cm  +  k  +  le 


■  y{n  -  m) 


*2(t./)-e  \x, (f,/)i2  y  i  c,c;c(,  +  *c;+*e-'(--") 

1  m—2  j 

„*-(*+ 1)  N-(k+D  .  .  ,  A 

+  £  !*,<8-t,/)|2  y  £ 

(1—1  (*l-l  J 

N-k  *-<*+ 1)  .  . 

+  E  XfirJ)  x;  (8  -  T ./)  y  £  cHcmmc,  +  " 

(i—i  m- 1 

N-U+l)  N-k  .  , 

+  £  Xf  (8  -  f ,/)  x;  ir,f)  y  £c.c;cB  +  *  +  ,c;  +  *e-^-"  +  I 

it— I  m—  I 


,  -  y(i i  -  m  +  I) 


-V 
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With  the  random  characteristics  of  an  aperiodic  code  sequence  previously  discussed,  it  is  easy  to  verify 
that  the  third  and  fourth  terms  in  the  right-hand  side  of  Eq.  (4)  are  zeros  since  the  coefficients  are 
mutually  ’exclusive”  at  all  times;  i.e.,  the  conditions  such  that  the  terms  c„c*c„  +  *c*  +  *  + ,  and 
c„c^cm  +  »  +  ic*  +  *  contain  no  odd  powers  of  the  parameters  for  1  <  /  <  N  -  k  do  not  exist.  For  the 
first  and  second  terms,  there  are  only  two  cases  resulting  in  nonzero  mean,  that  is,  (a)  k  -  0  and  (b) 
k  ;*0  and  n  —  m.  Therefore  Eq.  (4)  can  be  reduced  to 


X2  (r ,/)  -  < 


[\X,  £  e>m)  +  Ur,(8-f,/|J  *~T+V '*(*"•) 

)  *  - 1  *■  - 1  n  -  1 


-  0 


\X,{r,f)\2  -'"(*"•) 

«- 1 


+  Ur,<8  -  r,/)|2 


N  -  (k  +  1) 

£  e-^") 

a  —  1 


(5) 


;  it  *  0, 


or  after  simplification 


X2(rJ) 


l^tf./)l2  -(1(1~_g7)(l  -  'e"')  +  _  f,/)|2  (N  ~  1);* 

<  |jr,(r./)|2Ar2  +  |JT,<8  -  r./)|2  (N  -  1);*  -  0  and  /  -  PRF 
||jr,(r./)|2(Af  -  *)  +  |J T,(*  -  r,/)|2  (N-k-lhk  9*0. 


0  and  /  9*  PRF 


(6) 


In  the  above,  X,(t,f)  and  Xf(&  —  f ,/)  can  be  obtained  by  direct  integration  as  defined  in  Eq.  (3), 
yielding  |Sn(ir/(8  —  t)|  (8  -  f )  and  \Sa(wfi)\r.  Here  Sa(x)  —  sin  x/x as  commonly  defined.  As  a 
result,  the  ambiguity  function  of  a  binary-coded  signal  is  expressed  as  follows: 


XHrJ)  - 


Sb2(*/B) 

SaHvfi 8  -  t))(N  -  *)(8  -  f)2  +  Sa2(wfr)iN  -k  -  l)2;*  -  1.2 . N  -  1 


(7) 


where  r  —  *8  +  f ,  0  <r  <  8,  and  N  and  8  are  the  number  of  code  elements  and  the  pulse  width. 
Note  that  the  above  expressions  are  valid  under  the  Mind  speed  condition;  i.e.,  the  doppler  frequency  is 
the  same  as  the  PRF  or  an  integer  multiple  thereof. 


With  Eq.  (7)  normalized  to  f  —  f/8,  /  —  /V8,  and  X2(r,f)  —  X2(r,f)/N2B2,  the  responses  of 

F0.F1/j.F1.F3/3,  and  Fj,  doppler  filters  for  N  -  PRF  •  T  —  127  are  shown  in  Figs.  1  through  5.  The 

code  element  of  N  —  127  is  specifically  chosen  here  so  that  the  responses  can  be  compared  with  those 
obtained  in  Ref.  [S]  for  the  case  of  pseudorandom  waveforms  (see  Discussion  and  Conclusion).  To  be 
consistent  with  the  name  'ambiguity  function”  used  in  radar  systems,  the  term  "doppler  filter"  as 
observed  in  matched-filter-bank  implementation  (doppler-steered  pulse  compressor)  of  radar  detection 
is  adopted  here.  For  F,  filters,  the  doppler  shift  is  such  that  the  total  phase  shift  across  the  whole  code 
is  2w  L  In  the  special  case  of  zero  doppler  shift,  the  mainlobe  response  can  be  obtained  from  Eq.  (7) 
with  both  /and  k  set  to  zero.  That  is, 

X2(r,0)/ N2B2  -  (1  -  t)2  +  t2(N  -  1 )/ Af2;  0  <  f  <  1.  (8) 

The  response  reduces  to  half  of  its  peak  value  at  r  -  0.29  (or  r  -  0.298)  for  the  typical  case  of 

N  —  4800  (Fig.  6).  For  k  ^  0,  the  sidelobe  response  is: 

X2(r  ,0)/  N262  -  (1  -  r  )HN  -  k)/N 2  +  t2(N  -k  -  D/TV2;®  * \  ^ _  , 
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Fig.  1  —  Fq  doppler -filter  response  for  •  random  binary-phase 
code  with  the  code  element  N  -  127 


Fig.  2  —  Fyi  doppier-fllter  response  for  a  random  binary-phase 
code  with  the  code  element  N  -  127 
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Fig.  3  —  F|  doppier-fllter  responae  for  a  random  binary-phaae 
code  with  the  code  element  N  —  127 
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which  varies  approximately  from  36.8  to  39.2  dB  below  the  main  peak  for  the  first  eight  cycles.  Clearly, 
the  function  governed  by  Eqs.  (8)  and  (9)  converges  to  a  triangular  one  near  the  origin  as  N  becomes 
large.  This  is  the  trivial  case  for  an  infinitely  long  random  sequence.  In  the  special  case  of  zero  time 
delay  by  setting  *  —  0  and  f  -  0  in  Eq.  (7),  the  frequency  response  of  the  ambiguity  is  simply  a 
sinHw  fN6)/(ir/N8)2  function  (Fig.  7)  which  is  identical  to  the  spectrum  as  seen  in  pseudorandom 
waveforms  [6]. 


Fig.  7  —  Ambiguity  function  frequency  response  for  ■  random  binary-phase  code 
with  the  code  element  N  -  4800 


DISCUSSION  AND  CONCLUSION 

In  communications  and  radar,  the  emitted  signals  are  either  pseudorandom  with  a  very  long  cycle 
time  or  coded.  From  the  reception  viewpoint,  the  discretely  coded  waveform  to  be  cross-correlated  has 
a  finite  code  length  and  can  be  considered  aperiodic.  Without  a  priori  knowledge  of  the  emitted  signal 
in  applications,  one  may  consider  the  waveform  truly  random  and  replace  the  cross-correlation  process¬ 
ing  by  its  ensemble  average.  Herein,  from  a  statistical  sense,  the  ambiguity  function  of  a  random 
binary-phase-coded  signal  is  analytically  derived.  Although  the  signal  detection  relies  on  the  realization 
of  cross-correlation  processing,  the  limitation  of  the  transmitter  power,  the  threshold  of  the  received 
signal-to-noise  ratio,  and  other  practical  considerations,  the  characteristics  of  the  ambiguity  function 
derived  here  provide  predictive  measures  and  specification  requirements  for  practical  system  design. 

The  derived  ambiguity  function  in  an  analytical  form  can  also  be  used  as  a  baseline  in  relation  to 
other  deterministic/ pseudorandom  correlation  functions  of  interest.  For  example,  responses  of  F0,  F\, 
and  F2  doppler  filters  for  a  127-element  maximum-length  shift-register  code  were  investigated  in  Ref. 
[S]  (Fig.  8).  Indeed,  their  fundamental  properties  such  as  resolution,  processing  gain,  and  sidelobe 
level  can  be  easily  predicted  through  the  corresponding  doppler  filter  responses  of  a  random  binary- 
phase  code.  The  responses  of  various  doppler  filters  directly  obtained  from  Eq.  (7)  with  the  code  ele¬ 
ment  of  N  —  127  are  presented  in  Figs.  I  through  5.  ft  is  observed  that  the  responses  of  F\-  and  F2- 
filters  are  quite  comparable  to  those  obtained  by  transmitting  an  aperiodic  maximum-length  pseudoran¬ 
dom  sequence.  Note  that,  in  Fig.  8,  only  discrete  sample  points  were  approximately  calculated  and  the 
specific  code  sequence  was  generated  from  a  seven-stage  shift  register  with  an  initial  condition  of  21 7  in 
octal. 
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For  further  applications,  the  derivation  of  the  main  result  presented  here  may  include  the  case 
that  coded  sequences  have  more  complex  occurrence  statistics.  The  consideration  of  target  fluctuation 
which  causes  a  spread  in  time  and  frequency  can  also  be  included.  In  the  latter  case,  the  output  of  the 
matched  filter  is  then  the  emitter  scattering  function  convolved  with  the  ambiguity  function  derived 
here. 


Fig.  8  —  Responses  of  the  doppier  Alters  with  a  zero  doppler  target  for  a  127-element 
maximum-length  pseudorandom  binary  code 
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